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Abstract 

Theorems on zeroes of the truncated generating function in the 
complex plane are reviewed. When examined in the framework of 
a statistical model of high energy collisions based on the negative 
binomial (Pascal) multiplicity distribution, these results lead to maps 
of zeroes of the grand canonical partition function which allow to 
interpret in a novel way different classes of events in pp collisions at 
LHC cm. energies. 

1 Introduction 

Following previous work ,1, the multiplicity distribution (MD) generating 
function in the dummy variable 2;, 

00 

G{z) = Y,PnZ'' (1) 

n=0 

has been recognised to be related to the grand canonical partition function of 
a system of n particles in statistical mechanics when z is identified with the 
fugacity variable; P„ is here the probability of detecting n charged particles 
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in full phase space (X^^o-^™ ~ -'-) ^^'-^ linked to the canonical partition 
function. 

It has to be stressed that in a realistic experiment n will never become 
infinity. Let us call M the maximum finite value of n which can be detected 
in an experiment. The generating function G{z) in this case is reduced to a 
polynomial of degree M in z: 

M 

GM{z) = Y,PnZ''. (2) 

n=0 

Let us now consider the algebraic equation 

Gm{z) = (3) 

in the complex z-plane. In the first place, notice that none of the roots Zi 
(i = 1, . . . , M) of equation Q can be real and positive, since the P„ in © 
are positive for any n. In the second place, recall that the M-truncated 
generating function (as any polynomial) can be factorised in terms of its 
roots as follows: 

M 

GM{Z)=PMX{{Z-Zi). (4) 

i=l 

By applying these considerations to phase transitions of a lattice gas (and 
of the Ising model), C.N. Yang and T.D. Lee have found that the roots Zi 
lie on a circle centered at the origin of the complex 2;-plane. The circle is open 
in a small sector bisected by the positive real axis; no zero lies in that sector. 
When M is odd, at least one root is real: it will of course be on the negative 
side. Gm{,z) turns out to be an analytic function along the positive real axis 
for any M. A phase transition in the sense of a non analytic behaviour of 
the thermodynamical functions (like pressure for instance) can only occur at 
points of the positive real z axis which are "in the thermodynamical limit" 
accumulation points of the zeroes of the algebraic equation (jHl). For M — > cxd, 
one should expect in general that the zeroes become closer and closer to the 
real positive axis and that the distance between the zeroes vanishes. 

The first application of these ideas to particle production goes back to the 
70's jHI, but the field flourished again (results were obtained in the context 
a discrete approximation to QCD cascades ||4j) in 1995, when E. De Wolf 
j3] studied, using the JETSET Monte Carlo generator, the zeroes of the 
generating function of the n-charged particle multiplicity distribution in e^e~ 
annihilation in the y-rapidity interval \y\ < 0.5 at 1000 GeV cm. energy. The 
zeroes were found to lie on a circle of approximately unit radius centered at 
the origin of the complex z-plane. Intriguing questions were also raised on 
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the development of this approach in different models and in different classes 
of events. 

After recalling and extending the main theorems on the subject, the 
present paper explores the properties of the distribution of zeroes in high 
energy collisions in the framework of the weighted superposition mechanism 
of different classes of events: each class is described by a Negative Binomial 
(Pascal) MD, whose characteristic parameters are n, the average charged 
particle multiplicity, and k, linked to the dispersion D hj k = n^/{D'^ — n). 

2 Maps of zeroes in the complex plane and 
the class of power series distributions. 

In 1997, T.C. Brooks, K.L. Kowalsky and C.C. Taylor ^ applied the En- 
estrom-Kakeya theorem [Zj to the scaled M-truncated generating function 

M 

Gm(A^) = 5^c„^" (5) 

n=0 

with coefficients given by Cn = X"Pn and the parameter A defined as the 
minimum of the ratio of two consecutive c„ coefficients: 

A = r[n^i^] = mm[Pn/Pn+i] = Pn^J Pn^,^+i- (6) 

The coefficients c„, for n = 0, 1, . . . , M, are real and positive and satisfy 
the condition 

Co > ci > ■ ■ ■ > Cm 

required by the Enestrom-Kakeya (EK) theorem in order to ensure that 
Gm{^z) differs from zero in the complex 2;-plane for \z\ < 1. Equation 
(0) can be rewritten in terms of the variable w = 1/z as 

Gm(A^) = z'^'gMiw), (7) 

with 

M 

9M{w) = Y,w"'Cm (8) 

m=0 

and 

Cm = A*^ ™ PM-m- (9) 

By choosing now the maximum of the ratio between two consecutive Cm 
coefficients: 

A = r[nmax] = max[P„/P„+i] = 
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one finds that the coefficients satisfy once again the EK theorem: they 
are indeed real and positive and ordered following the rule Cq > ■ ■ ■ > > 
• • • > Cm- Therefore gM{w) differs from zero for \w\ < 1 and the zeroes of 
guiw) lie necessarily in the region \w\ > 1. 

Accordingly, all the zeroes of the truncated generating function Gm{z) 
fall in the annular region in the complex z-plane centered at the origin and 
limited by r[njnin] and r[nmax]- The positive real axis is of course excluded. 

The mentioned theorem turns out to be quite useful when applied to the 
class of power series distributions, defined as follows: 

Pn = CnVPo. (10) 

The corresponding M-truncated generating function is 

M M 

Gm{z) = AmPo ('nb^'z" = AmPo (11) 

n=0 n=0 

where Am = (Xlnlo-^™) ^ normalisation factor. In the polynomial (fTT|l . 
the term Pq can be neglected and bz replaced by the variable u, producing a 
new scaled generating function in u: 

M 



HM{u) = YcnU'' (12) 



n=0 



which — as we will see in the following — is of particular interest for classifying 
different classes of events described in terms of NB (Pascal) MD with different 
parameters. In fact, applying the above considerations to the generating 
function 

Gi™-(.).|:pr'--'-(^^^)' (13) 



n=0 

with 

nl[n + k) 

leads to the M-truncated version of G[^^](^): 



(14) 



M M n yj ^ ^ 

GTi^) = AT E = E n i^^i^) ^0^"^" (15) 

with b = Notice that in JTHI), G^^^^{z) has a singularity for z = 1/6. 
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By defining now in ()15p 



i=l 

and 



k + i 



hz = u, 

the new renormalised generating function, neglecting tlie inessential factor 
Pq, becomes 

Hir\u) = En (i + ^) = atw'u". (17) 

n=0 i=l ^ ^ n=0 

Notice that the not-truncated version H^^^^(u) of (fTTj) has a singularity at 
u = 1. Relation (fT7j) has been used in its not-truncated form in Ref. 
where c|f^^' has been identified with the n-particle canonical partition func- 
tion; H^^^\u) is in this framework the truncated grand canonical partition 
function. It follows that if parameter k is less than one (as in the case of 
the third class of events in pp collisions at 14 TeV cm. energy [H]), then the 
coefficients c\^^^ are all real and positive, and decrease starting from c\^^\ 
The conditions for the applicability of the EK theorem are therefore also here 
satisfied, and the zeroes of the scaled generating function (fTTj) lie all outside 
the circle of unit radius \u\ = 1. 

This situation should be contrasted with the application of the EK the- 
orem when the parameter k is greater than one (as in the case of the soft 
and semi-hard classes of events in pp collisions ISJ. In this case in fact, the 
coefficients c|f^^^ are in increasing order and the zeroes lie all inside the circle 
of unit radius. 

As indicated in [7j and shown in [Hj, for the (not-rescaled) NB (Pascal) 
MD Eq. (US)), one finds 



r\n\ 



i^^^' lfn + 1 



p[NB] b\n + k 



thus all zeroes of the algebraic equation of the M-truncated generating func- 
tion lie inside the annular region of the complex z-plane delimited by 

r[nmin] and rfn^ax], with 

. , If M -^-L\ r ^ J- J-/, 

r[TimaxJ = T TT—T < Pl < '"PminJ = TT = - < I) (19) 



1 /M+ 1\ , , ,,11,, 
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and 

Notice that r[n] — >■ ^ when n — >■ oo, and r[n] — 1 when k (i.e. 6 — 1). 
In the rescaled version ()17|1 . r[n] = it tends to 1 when ri — > oo with 

constant k. It is remarkable that for k = 1, r[n\ is also equal to one and all 
zeroes lie on the circle of unit radius, the u = \ point being excluded (see 
Figures 1 and 2). 

The role of parameter a in the case of the NB (Pascal) MD is underlined 
by the fact that the sum of the inverse of the zeroes of the algebraic equation 
Q is given by 



Recall in fact that 



i=l 

M 



Zi -TO 



1=1 

and also that the sum of (n — 1) by {n — 1) product of zeroes equals to 

Pm 

for example for M = 3 one gets 

ZiZ2 + Z1Z2, + ^22:3 = (-1)^^. 



It follows that 



^ {-if-'Pu-iM-D/PM Pi 



{-lf'P,/PM Po' 



Notice that for the NB (Pascal) MD, Eq. (j2H) becomes 

I , , nk , , 

- = -bk = = -a. 22 

Zi n + k ^ ^ 

i=l 

Equation (PT|) is indeed a general property of polynomials and is independent 
of the value of M. In the w-rescaled version, one has 

M 

E- = -^- (23) 
1=1 



In addition, it should be pointed out that the NB parameters of the coeffi- 
cients of the M-truncated generating function (fT3jl . namely um and /cm, can 
be calculated directly from the zeroes of equation 0. One finds: 



M 
i=l 



(24) 



and 



Z^i=l 1-zi J 
yM ' 



(25) 



In the limit M — »• oo, the parameters um and kM coincide of course with the 
standard parameters of the not-truncated distributions. 

One question still to be answered concerns the multiplicity of zeroes of 
the algebraic equation (jHl), both in general and in the specific case of the NB 
(Pascal) MD. Let us consider two polynomials S{z) and T{z) of degree m 
and n respectively: 



S{z) 



+ Sr 



and 



T{z) = toz"" + t^z""-' + --- + t. 



[so ^ 0) 



(to + 0) 



It is known that S(^z^ and Ti^z^ possess a not constant common factor iff the 
determinant 



R 



ST 



det 



So S\ 

So si 



Sm—\ "^m 
Sm—\ 











\ 












tl 

to 




to 



So si 

tn—l tn 
tn~l 



Sm—1 













\ 







to 



tn I 



is equal to zero. 

It follows that, called S' the ffist derivative of the polynomial S'(z) with 
respect to z^ S'(z) will have at least one double root iff Rs^s' = 0. In order 
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to be sure that the zeroes of a M-truncated generating function Gm{z) are 
all distinct, it is therefore sufficient to show that none of the M roots Zi is 
such that GM^Zi) = G' m{z,i) = 0. 

By rewriting Pi^^\n, k) for simplicity in terms of parameters a and b 
(see equation (fTT?|l for the definition of a), one has 

Pi^^](a,6) = ^a{a + b) . . . {a + bin - l))Po, (26) 
nl 

and by properly normalising Gi\j{z) one finds, by induction. 



P ANB^ a'^'lia + b)ia + 2&) . . . (a + Mb)f-'Po''-' 
RgT\g'T'-^^ • ^ ^ 

Since clearly R„[nb] ^,[nb] is always positive for a and b greater than zero, 

it follows that all the zeroes of the algebraic equation G^^^^\z) = are 
distinct. Since all zeroes of G^IJj^\z) = are distinct, M conditions are 
needed to uniquely determine the set of zeroes. For = 1 for instance, one 
should expect that 

MM M 
i=l i=l i=l 

In fact, from (fT^ 

M [NB] 
1=1 '-M 

The same limit is obtained for 



M / [NB] \ 2 [NB] 

' ~ \ JNB] ) ^ [NB] ' 
i=l \ ^M / ^M 

as well as for the other powers up to M. Notice that the same conditions 
are satisfied for M — cx) and constant k: consequently, also in this case the 
zeroes stay on the circle of unit radius, as will be discussed in the following. 

3 Maps of zeroes for different classes of events 
in the rescaled complex z-plane. 

We propose in this Section to apply the previously discussed theorems and 
remarks to the study, in the rescaled complex z-plane (i.e., in the complex u- 
plane), of maps of zeroes of the algebraic equations obtained by truncating at 
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integer value M the NB (Pascal) MD generating function. This distribution 
is used in a statistical model of multiparticle production to describe different 
classes of events (or substructures) in different collisions. These substructures 
are for instance 2- and 3-jets events in e^e~ annihilation or events without 
mini-jets (soft events) and with mini-jets (semi-hard events) in pp collisions. 
Recently P , in pp collisions, the occurrence of a third class of hard events at 
14 TeV cm. energy has been postulated also in the framework of the weighted 
superposition mechanism of different classes of events, each described by 
a NB (Pascal) MD with characteristic NB parameters. The occurrence of 
this third class of events was motivated by the fact that, by assuming both 
strong KNO scaling violations and QCD behaviour, extrapolations of NB 
parameters in the LHC energy domain implied, in the semi-hard sample of 
events, a surprising reduction in the average number of clans jH]. 

In order to deepen the analysis of the properties of collective variables 
in terms of the maps of complex zeroes and in order to study the possible 
occurrence of a phase transition, we focus our attention on the three men- 
tioned classes of events at 14 TeV cm. energy in pp collisions. Accordingly, 
the NB parameters used in the following are those described in kgoft = 

^ i^soft 40, fcsemj— ?tard 3.7, nisemi—hard 87, kfiaj-d 0. 1212, Ti/jurd 460. 

We find that for the soft and semi-hard components, the parameter k 
is greater than one and the zeroes stay inside the circle of unit radius in 
the complex w-plane (see Figures la and lb). For the hard component, the 
parameter k is smaller than one and all the zeroes lie outside the circle of 
unit radius |-u| = 1 (see Figure Ic). In order to make our statements more 
transparent, in Figures 2a, 2b, 2c the behaviour of the three classes of events 
in the region around the point u = 1 has been magnified for different cut- 
offs but fixed value of k. In Figure 3 the above maps are shown together 
close to the point u = 1, with the same fixed cut-off for three values of k 
corresponding to the three classes of events. 

Notice that for = 1 all the zeroes stay on the circle of unit radius with 
the u = 1 point excluded. In addition, when M ^ oo the distribution of 
the zeroes of the algebraic equation obtained from the truncated and rescaled 
generating function converges towards the u = 1 point for all classes of events 
with both k > 1 and A; < 1. 

As mentioned in the Introduction, in order to test the occurrence of a 
phase transition one has to find the behaviour of the zeroes of the grand 
canonical partition function Q in the thermodynamical limit, i.e., with the 
size of the system going to infinity. The discussion is based on the approach 
contained in [T], where the generating function G{z) has been related to Q 



by 



Q = [Cm 
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Relation (j^H|) implies that the zeroes of the M-truncated grand canonical 
partition function Qm can be obtained from the zeroes of the M-truncated 
generating function Gm simply by rescaling those zeroes by a factor b (which 
was identified in FT] — it should be stressed again — with the fugacity). 

The behaviour of the zeroes will then depend on the volume and the 
number of particles: when the MD is of NB type, as we are using now, 
it is sufficient to know the behaviour of the NB parameters. At present 
however, a realistic link between k and n and the system size (the volume) 
is under investigation, hence postponed to a forthcoming paper. Here we 
limit ourselves to explore what happens when the thermodynamical limit 
corresponds to M (and n) going to infinity with k fixed. In the first place 
we will check whether the EK theorem can be used to confine the zeroes to a 
narrowing region of the complex plane. Recalling Eq. (fTHj) . we notice that for 
k = 1 one has r[n] = |. Thus minimum and maximum radii coincide and all 
the zeroes of Gm{z) lie on a circle of radius 1/6, independently of the chosen 
cut-off M. Numerical studies show that the opening around the positive real 
axis is reduced when M — oo; this suggests that z = 1/6 is an accumulation 
point of zeroes. It means of course that in this scenario for k = 1, the point 
u = 1 is an accumulation point of zeroes of the grand canonical partition 
function: since u on the real axis coincides with the fugacity, only for 6 — > 1 
will a phase transition take place. For k > 1 the zeroes of the grand canonical 
partition function obey the following relation 

.,1,, .,M + 1 

^^FminJ = j:< \U\ < r[nmaxj = ^ _^ ^ (29) 

while for k < 1 the two limits are exchanged. Such a region does not shrink. 
We have to rely again on numerical root-finding, which still indicates the 
closing in of the zeroes to m = 1 (see Figures 1 and Figures 2), and at the 
same time shows that they are much closer together than the bound given by 
relation allows. This feature again points in the direction of the presence 
of a possible phase transition at 6 — > 1. 

In conclusion, the use of the rescaled complex variable u = bz marks 
visibly the different behaviour of the map of zeroes of the hard class of events 
with respect to the maps of the soft and semi-hard ones in our statistical 
model. In addition, it leads to identify the coefficients clf^'^' of Eq. (fTT)|l with 
the canonical partition function for a system of n particles, and Hm{u) with 
the truncated grand canonical partition function of the same system (in this 
framework Ci = k, i.e., the parameter k coincides with the canonical partition 
function for a system with one particle, see also P). It should also be pointed 
out that in the rescaled version of this approach the NB parameter k alone 
controls the geometry of the distribution of the zeroes. This is indeed a quite 
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general property which is expected to be vahd also in other collisions and 
especially in those collisions involving events with large numbers of particles, 
like heavy ion collisions. Finally, it has been shown that the NB parameter 
k is equal to minus the sum of the inverse of the zeroes, independently of the 
value of the cut-off M. 

All these facts suggest to investigate further the thermodynamical content 
of the NB parameter k, and especially its dependence on temperature and 
volume, in view of a possible phase transition in the 6 — > 1 limit. Being 
indeed the fugacity b = e'*/'^-^^ (with fi the chemical potential), one finds 
that fi goes to in the limit 6-^1: this is a situation which is reached in a 
high temperature and high density real boson gas and which is expected to 
occur in a quark-gluon plasma. 

References 

[1] A.Giovannini, S.Lupia and R.Ugoccioni, Phys.Rev.D 65, 094028 (2002) 

[2] C.N.Yang and T.D.Lee, Phys.Rev. 87, 404 (1952); 87, 410 (1952) 

[3] K.J. Biebl and J. Wolf, Nucl.Phys. B 44, (1972) 301 

[4] B. Andersson, G. Gustafson and J.Samuelsson, Nucl.Phys .B 463, (1996) 
217 

[5] E. De Wolf, in Proceedings of the XIV Int. Symp. on Multiparticle Dy- 
namics, Eds. A.Giovannini, S.Lupia, R.Ugoccioni, World Scientific, Sin- 
gapore 1995, p. 15 

[6] T.C.Brooks, K.L.Kowalski and C.C.Taylor, Phys.Rev.D 56:5857-5861, 
1997 

[7] M.Marden, Geometry of polynomials (Mathematical Surveys Number 3, 
American Mathematical Society, Providence, Rhode Island, 1966, p.l36- 
137) 

[8] A.Giovannini and R.Ugoccioni, Phys.Rev.D 59, 094020 (1999) and 
Phys.Rev.D 60, 074077 (1999) 

[9] A.Giovannini and R.Ugoccioni, Phys.Rev.D 68, 034009 (2003) 



11 



Im [u] 





, 0..75* 


• • . , ^ 




0.5 






0.25 


• 


T 




• 


0!.15 -0.5 


-0.25 


0.25 0.5 0.75 




-0.25 






-0.5 






' -S..75 , 





Fig. la: rescaled map of zeroes, in the complex w-plane, of the NB (Pas- 
cal) MD at 14 TeV cm. energy with ksoft = 7 and cut-off M = 80 (M has 
been taken equal to 2nsoft)- Zeroes lie inside the circle of unit radius \u\ — 1. 
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Fig. lb: rescaled map of zeroes, in the complex w-plane, of the NB (Pascal) 
MD at 14 TeV cm. energy with ksemi=hard = 3.7 and cut-off M = 174 (M has 
been taken equal to 2nsemi=/iard)- Zeroes lie inside the circle of unit radius 
\u\ = 1. 
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Fig. Ic: rescaled map of zeroes, in the complex w-plane, of NB (Pascal) MD 
at 14 TeV cm. energy with khard = 0.1212 and cut-off M = 600. Zeroes lie 
outside the circle of unit radius \u\ = 1. 
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Fig. 2a: the rescaled map of zeroes of the NB (Pascal) MD of Fig. la has 
been magnified around the point w = 1 for different cut-offs M (white and 
black triangles correspond to M = 80 and M = 320 respectively, and white 
and black squares to M = 160 and M = 640). As M increases, the point 
u — 1 is shown to become an accumulation point of zeroes from the inside. 
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Fig. 2b: the rescaled map of zeroes of the NB (Pascal) MD of Fig. lb has 
been magnified around the point u — 1 ior different cut-offs M (white trian- 
gles correspond to M = 174, white squares to M = 348 and black triangles 
to M = 696). As M increases, the point m = 1 is shown to become an accu- 
mulation point of zeroes from the inside. 
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Fig. 2c: the rescaled map of zeroes of the NB (Pascal) MD of Fig. Ic has been 
magnified around the point ti = 1 for different cut-offs M (white triangles 

and squares correspond to M = 600 and M = 700 respectively, and black 
triangles and squares to M = 800 and M = 920). As M increases, the point 
M = 1 is shown to become an accumulation point of zeroes from the outside. 
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Fig. 3: maps of zeroes of the NB (Pascal) MD are shown for different values 
of the parameter k and the same fixed cut-off M = 800. White triangles 
correspond to kgoft of Fig. la, white squares to ksemi-hard of Fig- lb and 
black triangles to khard of Fig. Ic. The full line describes the circle of unit 
radius \u\ — 1. The figure illustrates that — as discussed in the text — the 
point M = 1 is an accumulation point of zeroes from the inside (for A; > 1) 
and from the outside (for A; < 1.) 
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